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Abstract 

We present a four-neutrino model with three active neutrinos and one sterile neu- 
trino which naturally has maximal z/^ v-r oscillations of atmospheric neutrinos and 
can explain the solar neutrino and LSND results. The model predicts V(. — > Uj- and 
J^e t^fj, oscillations in long-baseline experiments with L/E ^ 1 km/GeV with ampli- 
tudes that are determined by the LSND oscillation amplitude and argument controlled 
by the atmospheric 6m'^. 



There is growing experimental evidence that neutrinos oscillate [|l|]. The long-standing 
solar neutrino deficit 0, 0, the atmospheric neutrino anomaly and the recent results 

from the LSND experiment on neutrinos from /i+ and tt^ decay can all be understood in 
terms of oscillations between two neutrino species. The challenge is to describe all oscillation 
phenomena within a single model, since resonant oscillations for the sun, oscillations for the 
atmosphere, and the LSND data each require a different neutrino mass-squared difference 
(5m^ to properly describe all features of the data P]. For example, if the atmospheric Sm^ scale 
is raised to the LSND scale, one forfeits the recently reported zenith-angle dependence of the 
atmospheric neutrino flux Alternatively, if the solar 6m'^ is raised to the atmospheric 6m'^ 
scale, one flnds that: (i) the reduction in the solar neutrino flux is energy-independent 0, 
and (ii) near-maximal Ug — or — mixing is necessary to describe the observed solar z/g 
suppression; but in the context of a 3 x 3 unitary matrix, such large mixing is inconsistent 
with the near maximal — mixing deduced from the atmospheric data. Then also 
if bw? > 10~^ eV^, large z/g mixing with any neutrino species is excluded by the recent 



CHOOZ reactor data [|TOl. Hence the large suppression of the solar flux can only result from 
resonance-enhancement, which requires a very small mass scale compared to those indicated 
by the atmospheric and LSND data. Since three-neutrino models can have at most two 
independent mass-squared differences Sm^, apparently not all the data can be explained 
with just z/g, z/^, and z/,-. 



A viable solution is to postulate one or more additional species of sterile light neutrino |]Tl 
without Standard Model gauge interactions (to be consistent with LEP measurements of 
Z vv [0) thereby introducing another independent mass scale to the theory. The latter 
approach has been used with some success in the literature W\ M. The constraints of 
big-bang nucleosynthesis give the constraint 

brn^A < IQ-'^ eV^ , (1) 

on the mass-squared difference and oscillation amplitude A = sin^ 29 for oscillations 

between a sterile neutrino and an active neutrino flavor ||15|| . 

In this Letter, we examine a four-neutrino model (three active plus one sterileQ ) which 

naturally has maximal z/^ — > z/^ oscillations of atmospheric neutrinos and which can also 
^ In principle, more than one sterile neutrino species can exist. However, only the particular linear 



explain the solar neutrino and LSND results. We begin with a brief discussion of the three 
classes of experiments and the neutrino mass and mixing parameters needed to explain them. 
We then present a mass matrix whose eigenvalues consist of a nearly degenerate neutrino 
pair at ~ 1 eV and a nearly degenerate pair at low mass, as illustrated in Fig. 1. We show 
how the existing data almost uniquely fixes the model parameters and strictly determines 
what new phenomenology the model predicts. We find that the new observable signature 
for the model (in addition to the oscillations already indicated by the data) is z/g ^ and 
z/g ^ oscillations for L/E ^ 1 km/GeV. We discuss the possibility that such signals can 
be observed in long-baseline neutrino experiments such as those using intense muon sources 



at Fermilab ^ or KEK and detectors at SOUDAN, GRAN SASSO, or AMANDA as the 
target. We find that the SOUDAN and GRAN SASSO possibilities would probe some of the 
possible Uf. —>■ oscillation region. We also compare this model to another four-neutrino 
mass matrix parameterization [|1^ which has been proposed to explain the data, and discuss 
their similarities and differences. 

LSND. The LSND experiment searches for — * z/g oscillations from /i"*" decay at rest 
(DAR) and for —>■ Vf, oscillations from vr"'" decay in flight (DIF). The DAR data has higher 
statistics, but the allowed regions for the two processes are in good agreement and suggest 

—>■ vacuum oscillation parameters that lie along the line segment described by 

0.030 eV^ 



0.3 eV^ < 5m^sND = 7^ ^ < 2.0 eV^ . (2) 

Larger values for 5ml^^^ are excluded at 90% C.L. by the BNL E-776 [0 and KARMEN ^ 
oscillation search experiments, and smaller values for 5mLgND excluded by the Bugey 
reactor experiment, which looks for v^, disappearance |jT9|. If only 99% C.L. exclusion is 
required, ^^«lsnd high as 10 eV^ is allowed for Alsnd — -0025; values of ^"^LSND > 10 
are excluded by the NOMAD experiment |2^, while values above 3 eV^ are disfavored by 
the r-process mechanism of heavy element nucleosynthesis in supernovae . 

Atmospheric. The atmospheric neutrino experiments measure and Vf. (and their an- 
tineutrinos) created when cosmic rays interact with the Earth's atmosphere. One expects 
about twice as many muon neutrinos as electron neutrinos from the resulting cascade of 



combination of sterile neutrinos that mixes with Vg, is phenomenologically interesting. 
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pion and other meson decays. Several experiments ^ ^ obtain a v^/ z/g ratio that is about 
0.6 of the value expected from detailed theoretical calculations of the flux |2^. The Super- 



Kamiokande experiment has collected the most data and a preliminary analysis indicates 
that their results can be explained as —>■ oscillations with |l|, ^ 

3 X 10"^ eV^ < 5ml^^^ < 7 x 10"^ eV^ , 0.8 < Aatm. < 1-0 , (3) 

with the high end of each range favored. Although Vg oscillations (where Vg is sterile) 

could in principle explain the atmospheric data, big-bang nucleosynthesis excludes this possi- 
bility unless the chemical potential of the neutrinos is modified Independent of flux 
normalization considerations, the z/^ oscillation channel is strongly disfavored by the 

zenith angle distributions of the data [Q. The recent CHOOZ z/g disappearance experiment 
also excludes Ue i^^i oscillations with large mixing A ^ 0.2 for dm^ > 10^'^eV^ [|TD|. 

Solar. The solar neutrino experiments measure z/g created in the sun. There are 
three types of experiments, z/g capture in CI in the Homestake mine, z/g — e scattering at 
Kamiokande and Super-Kamiokande, and z/g capture in Ga at SAGE and GALLEX; each is 
sensitive to different ranges of the solar neutrino spectrum and measures a suppression from 
the expectations of the standard solar model (SSM)0. Matter-enhanced MSW oscillations 



of z/g i^ii,i^T 0, or Us 1 25] are sufficient to explain the data. For z/g — > Us the allowed 



parameter region ||26[ is bounded by 



3.5 X 10^^ eV^ < 5m1^^ < 7.5 x 10~^ eV^ , 2.5 x 10"^ < A^^i. < 1.6 x 10"^ . (4) 

The allowed regions for small-angle MSW z/g — > z/^, z/,- solutions are slightly smaller. However, 
if the LSND data is to be explained by z/^ —>■ z/g oscillations at the relatively large mass scale 
indicated in Eq. |^, and the atmospheric data by z/^ — > z/,- oscillations with the much smaller 
scale of Eq. 0, then the solar neutrino data would appear to suggest z/g Ug since the 5m1^ 
scale, which is also given by Eq. ^, is not consistent with Eq. ^. It is necessary that the 
eigenmass mo associated predominantly with be heavier than the mass mi associated 
predominantly with z/g so that it is z/g rather than z/g that is resonant in the sun, which 
requires 

5mQi = mg — m\ > . (5) 
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Complete description of the data. In order to explain all of the above data, one needs 
a model which includes the three different two-neutrino solutions described above. The 
appropriate mass scales Sttiqi for MSW solar, 6m'^2 atmospheric, and Sm^i — Sm^^ ~ 
Sm'^i ~ Sm'^i for LSND oscillations are provided by the mass hierarchy ml ^ itlq <^ — m|. 
Taking the small-angle solar MSW solution, the required oscillation amplitude hierarchy is 

Mass matrix ansatz. To describe the above oscillation phenomena, we consider the neu- 
trino mass matrix 

^ el e?e2 ^ 



M = m 



ete2 es 



e| 



(6) 



V 63 1 el y 

presented in the (z/^, u^, z/^, i^r) basis. The mass matrix M contains five parameters (m, ei, 
£2, es, £4), just enough to incorporate the required three mass differences and two small oscil- 
lation amplitudes A^^ and A^s- The large amplitude A^,- does not require a sixth parameter 
in our model, because the structure of the v^-Vt submatrix naturally gives maximal mixing 
here (more on this below). We note that changing the position of €3 from the Mer element to 
the Me^ element would cause the z/^ to oscillate into Ve instead of into v^. If nonzero terms 
are introduced at both the M^t and M^.^ positions, then the physics changes: both and Vr 
would mix with v^, at the LSND scale, and the Ve ~ mixing angle is also affected^. Here we 
choose to take the minimal M needed to describe the data and determine the consequences. 

For simplicity, we have taken the mass matrix to be real and symmetric; then M is 
diagonalized by an orthogonal matrix U . Since U is real, there is no CP violation (which 
should be small anyway since observable CP-violation requires more than one large mixing 
angle, while the data seems to indicate just one). The are assumed to be small and of the 
same order of magnitude; phenomenologically they turn out to be within a factor of two of 
each other. 



^Thc other zero terms could be taken as nonzero without changing the phenomenology discussed here 
as long as they are small compared to e^. Inclusion of a small nonzero M^e term merely increases the 
tiny eigenmass mi, while a small nonzero Msp, or A/^r gives the sterile neutrino a larger but nevertheless 
unobservable mixing with and Vt. 
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The smallness of Me^/lM^g — Mgel = 62 is designed to yield the small-angle MSW solution 
for the sun. Simple changes in the 2x2 — z/g sub-block of M would allow us to also consider 
the large-angle MSW solar solution, but since large mixing of sterile with active neutrinos 
is disfavored by the solar data ||26| and big-bang nucleosynthesis [|15] for these values, 
we do not pursue this option here. 

To a good approximation, the eigenvalues of the mass matrix in Eq. 1^ are 



mo 



me 



1 5 



mi 



777,2 



3 = (1 T e| + |e^) 



(7) 



which shows the desired hierarchy. The small relative mass splitting of the heavier masses 
7772, ^3 is governed entirely by the parameter e|. Defining 5mf 



mj—m'j, the LSND 



oscillations are driven by the scale 7?7^ 



6m: 



21 



6m 



31 



6m: 



20 



57?7|q, the atmospheric 



Ur oscillations are determined by 6m 



32 



4777^64, and the solar 



z/o oscillations are 



determined by 6mli ~ m'^ef. The charged-current eigenstates are approximately related to 
the mass eigenstates by 



u 



\ ) 



I 



1 

£2 
-6263 



V e2e3( 



-^2 
1 

e3e4 



1 ^ 
1 

\_ 



1 ^ 
73^3 

1 

1 



/ V ^3 / 

Note that z/q and z^i couple 



Unitarity holds to first order in the ef UW = 1 + 0{e^ 
predominantly to z/^ and u^, respectively, as desired. The near-degenerate 1^2 and z/3 are 
seen to consist primarily of nearly equal mixtures of z/^ and u^-. These results are shown 
schematically in Fig. 1. 

Oscillation probabilities. With real-valued U, the vacuum oscillation probabilities are, in 



general, given by 



(5a/3 - 4 ^ UaiUiSiUajU/Sj siu^ A 



(9) 



where Aji = 6m'j-L/AE = 1.27{6m'^ /eV'^){L/km)/{E /GeV). For the mixing in Eq. I, the 
off-diagonal vacuum oscillation probabilities, to leading order in ej for each Aij and ignoring 
oscillations smaller than O^e'j), are given by 



(2 sin^ A21 + 2 sin^ A31 - sin^ A32) 



+ e^e^2 sin^ Aso + 2 sin^ Ago - 4 sin^ Aoi) , (10) 

P{ue^Ur) - sin^ A32 + 2e^e2 (^sin^ A21 - sin^ A31) , (11) 

P{u^^Ur) - sin^ A32 + 2e2e2 (sin^ A31 - sin^ A21) , (12) 

P(z/e^z/,) ~ 4e2sin2Aoi, (13) 

P(z/,. ^z/,) - Aejelsin^Ao^, (14) 

where Aqi ^ A32 ^ A20 — A30 — A21 — A31 due to the spectrum of the neutrino mass 
eigenvalues. 

For small L/E only the leading oscillations A20 — A21 — A30 — A31 contribute, and the 
only appreciable oscillation probability is 

P(z/, -^z/^)~4e2sin2A, (15) 

where A = m?L/4:E. From Eq. |15| we can fix two model parameters 

<^"^LSND = ) ^LSND = 4e3 . (16) 

The vacuum oscillation length associated with the LSND 5m?' scale is 

= 47rE/W = 2.5 km(E/GeV)((5m^sND/eV2)-^ . (17) 

For L/E typical to atmospheric or long baseline neutrino experiments, the oscillations in 
A assume their average values. The A32 oscillation is now evident, and the non-negligible 
oscillation probabilities in vacuum are 

^32) , (18) 

(19) 
(20) 

From Eq. ^ 

^^Itm. = ^"^32 - 4m2e2 , Aatm. = 1 , (21) 

which determines another parameter of the model. The model automatically gives maximal 
i^fj. — ^ oscillations for atmospheric neutrinos, while oscillations in other channels are sup- 
pressed. The v^-Vt maximal mixing is natural in the sense that it results from the large 
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P{ 






~ £3 (2 - sin 


Pi 




Vr) 


~ el sin^ A32 


P{ 


- 




~ sin^ A32 . 



value of the M^^ matrix element relative to the diagonal M^^ and M^-^ elements, without 
any need for fine tuning of the difference |M^^ — M^,-|. The oscillation length resulting from 
the 5m|2 scale is 

= 500 km(E/GeV)(5m^t^/5 x 10" W)"^ . (22) 

Finally, for very large LjE ^ ('^"^atm./^^^)^^ km/GeV, sin^ A32 averages to | and the 
appreciable oscillations in vacuum are (to leading order in the e^) 









Ae\ sin^ Aqi , 


(23) 


P{^e- 








(24) 








2 3 ' 


(25) 




Z^r) 




1 

2 ' 


(26) 


P{^,- 






Aelel sin^ Aqi . 


(27) 



The solar data can then be explained with the usual MSW matter-enhanced mechanism 
(including the proper sign of ^mgi in Eq. ^ if the parameters in vacuum satisfy 

5m2„, = 5ml^ ~ Amh\ , A,,,, = Ael . (28) 

Summarizing the above analysis, the model parameters are related to the observables by 

For the specific values Sml^^Y) = 2 eV^, Alsnd = 2.5 x 10~^, Smlf.^ = 5 x 10"^ eV^, 
^atm. = 1, (^""^soi. = 4 X 10~^ eV^, and Agoi = 1 x 10"^, we obtain 

m = 1.4eV, 61 = 0.038, €2 = 0.050, es = 0.025 , €4 = 0.025. (30) 

The corresponding neutrino mass eigenvalues are (in eV) 

1 1 
mo = 2xlO"^ ml=4xlO"^ -(ma + ma) ~ 1.4 , -(ms - ma) ^ 9 x 10"^ . (31) 

For these masses J2 ~ 3 eV, which according to recent work on early universe formation 
of the largest structures provides an ideal hot dark matter component [^ . 
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If instead we use the lowest allowed mass scale for the LSND experiment we obtain 
LSND — 0-3 eV^ and Alsnd = 4 x 10~^, in which case 

m = 0.55eV, ei = 0.060 , es = 0.050 , 63 = 0.10, €4 = 0.065, (32) 



with corresponding mass eigenvalues (in eV) 

mo = 2xlO"^ ml = 2xl0~^ ^(mg + ma) ~ 0.55 , ^(mg - ma) ~ 2.3 x 10"^ . (33) 

In either of the above examples, the 5m^ scale for the atmospheric neutrino oscillation can 
be adjusted simply by varying e^. Also in either case, the two heaviest masses provide rehc 
neutrino targets for a mechanism that may generate the cosmic ray air showers observed 
above > lO^o eV 

Model predictions. The model is constructed to provide the effective two-neutrino oscil- 
lation solutions for the LSND, atmospheric and solar data. The Solar Neutrino Observatory 
(SNO) ||31|, which can measure both charge-current (CC) and neutral-current (NC) inter- 
actions, will be able to test the Ue — > i^s solar oscillation hypothesis: in the sterile case the 
CC/NC ratio in SNO would be unity and both CC and NC rates would be suppressed from 
the SSM predictions. 

Given the order of magnitude of the Sm^j and Uaj, observable new phenomenology occurs 
for L/E ^ 1 km/GeV in the oscillation channels 

P{l^e Z^m) - ^^LSNd(2 - Sin^ Aatm.) , (34) 

1 



P(z/e l^r) - ^^LSND siu Aatm. , (35) 

where Alsnd ~ 0{1%) is the oscillation amplitude which describes the LSND results and 
Aatm. = l.27Sml,^ L/E ~ (^m^^^ /5 x lO-^ eV^) (L/ 157km) {GeV / E) is the oscillation argu- 
ment which describes the atmospheric neutrino data. In addition to the z/^ Ve oscillations 
due to A in Eq. |T^, which reach their oscillation-averaged value of |v4lsnd, the model pre- 



dicts new oscillations in the ^fi and z/g — >■ 1/^ channels with common oscillation length 
determined by Aatm. and amplitude given by |Alsnd- 



How can the oscillation probabilities in Eqs. ^ and ^ be tested? A list of experiments 
currently underway or being planned to test neutrino oscillation hypotheses is given in Ta- 
ble 1 [0. In each case the oscillation channel and the parameters which are expected to be 
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tested are shown. Many of these experiments will not provide any constraints on the new 
phenomenology, although many provide some check on the existing LSND or atmospheric 
neutrino results (those that provide tests are noted in the table). The KARMEN upgrade, 
Booster Neutrino Experiment (BooNE at Fermilab), ORLANDO at Oak Ridge, and MINOS 
(Fermilab to SOUDAN) can test the LSND —>■ z/g oscillations, and for the experiments 
that probe significantly below 1 eV^, may be able to detect the contribution of the 
additional oscillation due to the sin^ A32 term in Eq. MINOS and ICARUS also aim to 
detect Ur and should be able to probe some of the atmospheric neutrino allowed region for 
Ur. NOMAD, CHORUS, and TOSCA at CERN and COSMOS at Fermilab will test 
Ut- oscillations at most down to ^0.1 eV^ for maximum amplitude; these do not 
probe our model as there are no appreciable z/^ —>■ Vr oscillations in that region. Reactor 
experiments at Palo Verde in Arizona and with the BOREXINO detector in Europe will test 
Ve disappearance involving appreciable mixing angles, but will not test our model since the 
largest Ve vacuum oscillations are the A = .04 level or less. 

Long-baseline experiments with an intense Ue or Ug neutrino beam which can detect r's, 
and hence see the —>■ Vt oscillations in Eq. |3^, can provide a definitive test of the new 
phenomenology of our model. High intensity muon sources |T6[ can provide simultaneous 
high intensity and Ve (or and Ve for antimuons) beams with well-determined fluxes, 
which could then be aimed at a neutrino detector at a distant site. It is expected that 
r's will be detected through their /i decay mode and that a charge determination can be 
made, so that one can tell if the r originated from — > or z/g oscillations. Current 
proposals [|l|] consider SOUDAN (L = 732 km) or GRAN SASSO (L = 9900 km) as the far 
site from an intense muon source at Fermilab. These experiments could also observe 
oscillations via detection of "wrong-sign" muons. The neutrino energies are in the 10-50 GeV 
range. Assuming that low backgrounds can be achieved, the sensitivity to is roughly 
proportional to the inverse square root of the detector size (given the same neutrino energy 
spectrum at the source); the 5rn? sensitivity does not depend on detector distance L because 
although the flux in the detector falls off with L^, the oscillation argument grows with for 
small 5m?L/ E. For 20 GeV muons at Fermilab and a 10 kT detector at either SOUDAN or 
GRAN SASSO, the single-event sensitivity for Uj. oscillations is about 8 x 10^^ eV^ 



10 



for maximal mixing [16|. For large Sm^, the oscillation amplitude single-event sensitivity is 



roughly inversely proportional to the neutrino flux at the detector divided by the detector 
size; about 6 x 10^^ for SOUDAN and lO^^ for GRAN SASSO pf. In general, the closer 
detector has comparable Sm"^ sensitivity but better A sensitivity. 

Our model predicts z/g —>■ Vr oscillations with amplitude Alsnd/4 (which ranges from 
0.0025 to 0.04) and mass-squared difference of 5m1^^ (which ranges from 3 x 10~^ to 7 x 10^'^ 
eV^). The region of possible Vf. Uj. oscillations in our model and the regions which can be 
tested at the SOUDAN and GRAN SASSO sites are shown schematically in Fig. 2, along with 
the favored parameters for the LSND, atmospheric neutrino, and solar neutrino oscillations. 
Such experiments would be sensitive to some of the z/g — > region, though they may not 
cover the low-mass, small-amplitude part. These searches would also be able to test the 
oscillations in Eq. and the atmospheric z/^ — ^ Vr oscillations. Additionally, long 
baseline experiments to the AMANDA |Q detector from Fermilab or KEK may be useful 
in probing oscillations with small 5m?. 

Neutrinoless douhle-fi decay. From the form of U and the mass eigenvalues one can readily 
see that neutrinoless double-/? decay is unobservable in our model. If neutrinos are Majorana 
particles, then the {PP)ou decay rate is proportional to 



<m^>=\Y, Uejrrijl ~ mel < 10~^eV , (36) 
j 

which is well below the present limit of ~ 0.5 eV P^, and less than improved bounds 



realizable in the future. Note that possible CP-violating relative Majorana phases which we 
have ignored in our model can give smaller < m^, > via a cancellation in the leading terms, 
but cannot give larger < >. 

Hot dark matter. The contribution of the neutrinos to the mass density of the universe 
is given by = J^f^u/ih'^Q^ eV), where h is the Hubble expansion parameter in units of 
100 km/s/Mpc p3; with h = 0.65 our model implies Q^, ~ 0.05. An interesting test of 



neutrino masses is the Sloan Digital Sky Survey (SDSS) [^. For two nearly degenerate 
massive neutrino species, sensitivity down to about 0.2 to 0.9 eV (depending on Q and h) is 
expected, providing coverage of all or part of the LSND allowed range (m =0.55 to 1.4 eV 
in our model). 
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Resonant enhancement in matter. The curves in Fig. 2 assume vacuum oscillations. In 
general, large corrections to oscillations involving v^. and Vs are possible due to matter; 
the Ve diagonal element in the effective mass-squared matrix receives an additional term 
2\/2GfN(,E from the CC interaction, and the Ug diagonal element receives the contribution 
\/2GFNnE (relative to the active neutrinos) because it does not have NC interactions, where 
Ne and Nn are the electron and neutron number density, respectively. In our model, however, 
these corrections do not significantly affect the large m\ and m\ mass eigenvalues as long as 
E <^ 1 TeV, and hence only modify the ^wiqi oscillation argument. For E > Q.l GeV, the 
only significant change in the mixing parameters is that the z/g — Vs mixing 62 in Eqs. |lO| to 
p!4| is suppressed. The result is that for small and intermediate L/ E {i.e., all experiments 

, |TB|, and ^Op, matter does not appreciably change the observable 



described by Eqs. 



phenomenology of the model. For large L/ E, such as when E' ^ 10 MeV for solar neutrinos, 
there can be maximal z/g — Vs mixing which changes Eqs. |23, |23, and to 



sin^ A 



01 



sin^ A 



01 



el sin^ Aoi 



(37) 
(38) 
(39) 



In this case the only significant effect of matter (other than the MSW enhancement of 
i^e that leads to the solar neutrino suppression) is to enhance the Ug oscillations 

and introduce a new oscillation in the Ue —>■ z/^ channel, although the amplitude of these new 
oscillations never gets above about 10^^. Hence we conclude that the matter corrections for 
the mass matrix in Eq. ^ probably have no observable consequences. 

Other models. Are other viable neutrino mixing schemes possible? A different form for 
the neutrino mass matrix is 



M = m 



el e?62 



V 



e?e2 





es 
e2 








(40) 



This alternate form contains one more parameter, ef, than the mass matrix in Eq. |^. Fine- 
tuning of this additional parameter is necessary to achieve maximal mixing in the 1^2 — V3, 
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sector. Eq. 50 is a generalization of the particular form used in Ref. [0 which has = 2e|. 
Again, zero elements can be taken nonzero as long as they are very small. The eigenvalues 
are given approximately by 



mo 



-1 ? 



nil 



-me 



3 ' 



"^2,3 



m 



+ ei 



(41) 



The charged current eigenstates are approximately related to the mass eigenstates by 
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(42) 



where 



and 



/?-e?/(e? + e^) 



s = sin 9nr = —F^ 



1 + 



1/2 



(43) 



(44) 



determines the z/^ — z/^ mixing and c = cos 6'^,-. We also have 

(^"^L. = - ~ m^{e\ - e^) , A.^;. ~ 4/3^6^ . 



(45) 



Note that m\ cannot be larger than mo, since this would lead to a negative ^m?^^\ and would 
imply that z/g and not z/g undergoes an MSW enhancement, contrary to the data. Hence ei 
must be larger than 63 with this matrix form, which in turn limits /5 to the interval [0.5, 1]. 
Then £2 must be raised by the factor relative to the matrix form in Eq. ^ to give the 
proper MSW mixing. The LSND constraints are the same as in Eq. |T6[ In the z/^ — z/,- sector, 
the parameters are determined by 



5m^,^.^2mV(ei-ei)2 + 46|, A, 



4et 



(46) 



- (ei - ei)2 + 4e| 

In this scenario the ea term significantly affects not only the lightest mass eigenvalue but also 
the mass splitting and mixing angle of the two heavy states. Consequently, some fine-tuning 
is necessary to achieve the proper phenomenology. Maximal mixing occurs only when ea ~ €5, 
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in which case Sn^l^^^ = Am?e\ just hke our previous scenario. However, in the absence of 
such fine-tuning, submaximal mixing is probable, in which case a different value for €4 is 
required to generate the correct 5m1^^ . 

We may again solve for the parameters directly in terms of the observables. In the present 
case there are six parameters and six observables including the (now) possibly non-maximal 
amplitude for atmospheric oscillations. The result is 



2 ^sol. 

^2 = —T 



e3 = ^LSND/4, (48) 



^^^5^^ ' 5~ 3 = ±2Vl-^atm.T--2 • (49) 



Since the matrix form in Eq. ^ requires some fine-tuning to explain the data, some higher 
order terms must be retained in the expressions for the parameters. 

Using the same input parameters as before, including maximal mixing in the atmospheric 
neutrino experiments (which implies €5 = €3), we find for the largest ^w-lsnd solution 

m = 1.4eV, 61 = 0.039, 62 = 0.070, 63 = 65 = 0.025, 64 = 0.025, (50) 

with mass eigenvalues (in eV) 

mo = 2.1xlO~^, mi=0.9xl0"^, i(m3 + 7/12) ~ 1.4 , ^(77^3 - ms) ~ 9 x 10^^ . (51) 

For the smallest ^m^gND solution, we obtain 

m = 0.55eV, 61 = 0.103, 62 = 0.070, 63 = 65 = 0.100, 64 = 0.065, (52) 

with masses (in eV) 

1 1 

mo = 5.8xlO"^ mi = 5.5xl0~^ -(mg+ms) ~ 0.55 , -(m3-m2) ^ 2.3x 10"^ . (53) 

In Eq. ^ some fine-tuning between 65 and 63 (to the 3% level) is needed for 5m1^^ to have 
the correct sign and magnitude. In either Eq. ^ or ^ the mass scale for the atmospheric 
neutrino oscillation can also be adjusted simply by varying 64 (for maximal mixing), or 
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by adjusting €4 and €5 (for non-maximal mixing). The only new phenomenology is again 
for LjE \ km/GeV, except that sin^ Aatm. in Eqs. ^ and ^ is now replaced 
by v4atm. sin^ Aatm. when Aatm. 7^ 1- The possible Vf. oscillation amplitude is reduced by 
a factor Aatm. (which is apparently 0.8 or higher), which shifts the predicted region in Fig. 2 
slightly to the left; otherwise this model is very similar to the model of Eq. ^ 

Summary. In this letter we have presented a four-neutrino model with three active neutri- 
nos and one sterile neutrino which naturally has maximal z/^ — » Uj. oscillations of atmospheric 
neutrinos and can also explain the solar neutrino and LSND results. The model predicts 
z/g — * Ur and Ug oscillations in long-baseline experiments with L/E ^ 1 km/GeV 

with amplitudes that are determined by the LSND oscillation amplitude and Sm^ scale de- 
termined by the oscillation scale of atmospheric neutrinos. Neutrino beams from an intense 
muon source at Fermilab or KEK with a detector at the SOUDAN or GRAN SASSO sites 
may be able to test part of the parameter region for these oscillations channels. 
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Table 1: Current and planned neutrino oscillation experiments. Check- marks denote acces- 
sible oscillation channels. The Sm^ and sin^ 29 sensitivies are given. 
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Figure 1: Neutrino mass spectrum, showing the approximate flavor content of each mass 
eigenstate, and showing which mass sphttings are responsible for the LSND, atmospheric, 
and solar oscillations. 
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Figure 2: Predicted region in the effective 5m^-sin^ 29 parameter space for z/g Ut oscilla- 
tions in the four-neutrino model (solid rectangle), which is determined by ^ of the LSND 
^ oscillation amplitude and the atmospheric neutrino —>■ oscillation Sm^ scale. 
The dashed curves show the potential limits that can be set by neutrino beams from an 
intense muon source at Fermilab [|I^ to detectors at the SOUDAN and GRAN SASSO sites 
for muons with energy of 20 GeV. Also shown are the parameters for the solar Ue — >■ z/^ 
oscillation. 
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